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Experimental study of pprj dynamics in the pp — > ppr] reaction* 
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A high statistics measurement of the pp — > ppr) reaction at an excess energy of Q = 15.5 MeV 
has been performed at the internal beam facility COSY-11. The stochastically cooled proton beam 
and the used detection system allowed to determine the momenta of the outgoing protons with a 
precision of 4 MeV/c (a) in the center-of-mass frame. The determination of the four- momentum 
vectors of both outgoing protons allowed to derive the complete kinematical information of the rypp- 
system. An unexpectedly large enhancement of the occupation density in the kinematical regions 
of low proton-77 relative momenta is observed. A description taking the proton-proton and the 
77-proton interaction into account and assuming an on-shell incoherent pairwise interaction among 
the produced particles fails to explain this strong effect. Its understanding will require a rigorous 
three-body approach to the pprj system and the precise determination of contributions from higher 
partial waves. We also present an invariant mass spectrum of the proton-proton system determined 
at Q = 4.5 MeV. Interestingly, the enhancement at large relative momenta between protons is visible 
also at such a small excess energy. 

In contrast to all other determined angular distributions, the orientation of the emission plane 
with respect to the beam direction is extracted to be anisotropic. 

PACS numbers: 13.60.Le; 13.75.-n; 13.85.Lg; 25.40.-h; 29.20.Dh 



MANIFESTATION OF THE 77-P-P 
INTERACTION 



Due to the short live time of the flavour-neutral 
mesons (eg. 7r°, 77, 77'), the study of their interaction with 
nucleons or with other mesons is at present not feasible in 
direct scattering experiments. One of the methods per- 
mitting such investigations is the production of a meson 
in the nucleon-nucleon interaction close to the kinemati- 
cal threshold or in kinematics regions where the outgoing 
particles possess small relative velocities. A mutual in- 
teraction among the outgoing particles manifests itself 
in the distributions of differential cross sections as well 
as in the magnitude and energy dependence of the total 
reaction rate. 

In the last decade major experimental 0, 0, H 0, IE IE 
HQ and theoretical 0, Ea E3, E3, efforts were concen- 
trated on the study of the creation of 7r^, 77 and 77' mesons 
via the hadronic interactions 0, 

EE 

lid . Measurements 
have been performed in the vicinity of the kinematical 
threshold where only a few partial waves in both initial 



and final state are expected to contribute to the produc- 
tion process. This simplifies significantly the interpreta- 
tion of the data, yet still appears to be challenging due 
to the three particle final state system with a complex 
hadronic potential. The determined energy dependences 
of the total cross section for <q' HQ and 77 0,11 II 110 
mesons in proton-proton collisions are presented in fig- 
ure n Comparing the data to the arbitrarily normal- 
ized phase-space integral (dashed lines) reveals that the 
proton-proton FSI enhances the total cross section by 
more than an order of magnitude for low excess energies. 
One recognizes also that in the case of the rf the data 
are described very well (solid line) assuming that the on- 
shell proton-proton amplitude exclusively determines the 
phase-space population. This indicates that the proton- 
77' interaction is too small to manifest itself in the exci- 
tation function within the presently achievable accuracy. 
In case of the 77 meson the increase of the total cross sec- 
tion for very low and very high energies is much larger 
than expected from the 1 5'o final state interaction be- 
tween protons (solid line), though for both the pp — ► pprj 
and pp — ► ppr/' reactions the dominance of the 3 Pq — ¥ Sqs 
transition 1 is expected up to an excess energy of about 
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The transition between angular momentum combinations of the 
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40 MeV and 100 MeV, respectively 14] . The excess at 
higher energies can be assigned to the significant onset of 
higher partial waves, and the influence of the attractive 
interaction between the rj meson and the proton could be 
a plausible explanation for the enhancement at thresh- 
old. A similar effect close to threshold is also observed in 
the photoproduction of 7/ via the 7c? — > pnr\ reaction [l9| 
indicating to some extent that the phenomenon is in- 
dependent of the production process but rather related 
to the interaction among the rj meson and nucleons in 
the Sn(1535) resonance reg ion. Indeed, a simple phe- 
nomenological treatment |l3tll4ll2"i| - based on factoriza- 
tion of the transition amplitude into the constant primary 
production and the on-shcll incoherent pairwise interac- 
tion among outgoing particles - describes very well the 
enhancement close to the threshold (dotted line). How- 
ever, this approach fails for the description of the invari- 
ant mass distribution of the proton-proton and proton-77 
subsystems determined recently at Q = 15 MeV by the 
COSY-TOF collaboration [Ij. The structure of this in- 
variant mass distributions, which we confirm in this pa- 
per utilizing a fully different experimental method, may 
indicate a non-negligible contribution from the P-waves 
in the outgoing proton-proton subsystem 1 1 1 1 ] . These can 
be produced for instance via 1 Sq — > 3 Pqs or 1 £>2 -^ 3 P2S 
transitions. This hypothesis encounters, however, diffi- 
culties in describing the excess energy dependence of the 
total cross section. The amount of the P-wave admixture 
derived from the proton-proton invariant mass distribu- 
tion leads to a good description of the excitation function 
at higher excess energies while at the same time it spoils 
significanly the agreement with the data at low values 
of Q, as depicted by the dashed-dotted line in figure ^ 
However, these difficulties in reproducing the observed 
energy dependence might be due to the particular model 
used in reference 0, and thus higher partial wave con- 
tributions cannot a priori be excluded. In contrast to the 
P-wave contribution the three-body treatment [l8| of the 
pprj system (dashed-double-dotted line) leads to an even 
larger enhancement of the cross section near threshold 
than that based on the Ansatz of the factorization of the 
proton- proton and proton- 77 interactions. It must be kept 
in mind however that a too strong FSI effect predicted 
by the three-body model must be partially assigned to 
the neglect of the Coulomb repulsion in this preliminary 



initial and final states are described according to the conventional 
notation :8j in the following way: 

where, superscript "i" indicates the initial state quantities. S 
denotes the total spin of the nuleons, and J stands for the overall 
angular momentum of the system. L and I denote the relative 
angular momentum of the nucleon-nucleon pair and of the meson 
relative to the NN system, respectively. The values of orbital an- 
gular momenta are commonly expressed using the spectroscopic 
notation (L = S,P,D,..., and 1 = s,p,d,...). 




Q [ MeV] 

FIG. 1: Total cross section for the reactions pp — > pprj' (cir- 
cles) and pp — » ppr\ (squares) as a function of the centre-of- 
mass excess energy Q. Data are from refs. Q, 0, H S ■ 
The dashed lines indicate a phase-space integral normalized 
arbitrarily. 

The solid lines show the phase-space distribution with in- 
clusion of the So proton-proton strong and Coulomb in- 
teractions. In case of the pp — > pprj reaction the solid line 
was fitted to the data in the excess energy range between 15 
and 40 MeV. Additional inclusion of the proton-?? interac- 
tion is indicated by the dotted line. The scattering length 
of a pv = 0.7 fm + i0.4fm and the effective range parame- 
ter r PV — — 1.50 fm — i0.24fm 17] have been arbitrarily cho- 
sen. The dashed-dotted line represents the energy dependence 
taking into account the contribution from the 3 Po — » 1 5os, 
1 5o —> 3 Pos and 1 D2 — > 3 P2S transitions [Tljl . Preliminary re- 
sults for the 3 Po —* 1 Sos transition with the full treatment of 
the three-body effects are shown as a dashed-double-dotted 
line [l q . The absolute scale of dashed-double-dotted line was 
arbitrary fitted to demonstrate the energy dependence only. 

calculations 0]. These illustrate that the simple phe- 
nomenological approach shown by the dotted line could 
fortuitously lead to the proper result, due to a mutual 
cancelation of the effects caused by the approximations 
assumed in calculations and the neglect of higher partial 
waves. This issue will be discussed further in section llVl 
after the presentation of the new COSY-11 data. The 
above considerations show unambiguously that for the 
complete understanding of the low energy pprj dynamics, 
in addition to the already established excitation func- 
tion of the total cross section, a determination of the 
differential observables is also necessary. These will help 
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FIG. 2: Definition of the centre-of-mass kinematical variables used in this article for the description of the ppjy system. In 
the centre-of-mass system, the momenta of ejectiles lie in an emission plane. Within this plane the relative movement of the 
particles is fixed by the square of the invariant masses s pp and s PV . As remaining three variables needed to define the system 
uniquely we use either <j>l, 6^, and ip shown in the upper panel or (j)* N , 9* N , and ipN defined in the lower panel. TV is a vector 
normal to the emission plane, which can be calculated as the vector product of the centre-of-mass momentum vectors of the 
outgoing protons. As an example two extreme orientations of the emission plane are shown in the right-lower panel. For further 
descriptions see the text. 



to disentangle effects caused by the proton-77 interaction 
and the contributions from higher partial waves. In this 
article we present distributions determined experimen- 
tally for two sets of orthogonal variables fully describing 
the pp?7 system, which was produced at an excess energy 
of Q = 15.5 MeV via the pp — > yyr\ reaction using the 
COSY-11 HHH facility at COSY [2$. 

II. CHOICE OF OBSERVABLES 

For the full description of the three particle system 
five independent variables are required. In the center- 
of-mass frame, due to the momentum conservation, the 
momentum vectors of the particles are lying in one plane 
often referred to as the emission, reaction or decay plane. 
In this plane (depicted in figure |2J) a relative movement 
of the particles can be described by two variables only. 
The square of the invariant masses of the di-proton and 
proton- r) system denoted as s pp and s pv , respectively, 



constitute a natural choice for the study of the interac- 
tion within the pp?7 system. This is because in the case of 
non-interacting objects the surface spanned by these vari- 
ables is homogenously populated. The interaction among 
the particles modifies that occupation density and in con- 
sequence facilitates an easy qualitative interpretation of 
the experimental results. 

The remaining three variables must define an absolute 
orientation of the emission plane in the distinguished co- 
ordinate system. This may be realized for example by 
defining the orientation for the momentum of the arbi- 
trarily chosen particle in the center-of-mass frame and 
the angle which describes the rotation around the direc- 
tion fixed by that particle. In one of our choices made 
following reference |25l |26| the corresponding variables 
are the polar and the azimuthal angle of the r\ momen- 
tum vector, depicted in figure[2]as (f>* and 6*, respectively, 
and the angle ip describing the rotation around the di- 
rection established by the momentum of the r\ meson. 
Figure [21 demonstrates that such rotation neither affects 
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the rj meson momentum nor changes the configuration of 
the momenta in the emission plane. In the case of ex- 
periments with unpolarized beams and targets the only 
favoured direction is that of the beam line. Therefore, as 
a zero value of the tj) angle we have chosen the projection 
of the beam direction on the plane perpendicular to the 
momentum vector of the rj meson. Note that we identify 
the z-axis with the beam direction. The </>*, 9* and 
variables can also be interpreted as Euler angles allowing 
for the rotation of the emission plane into a xz-plane. 

As a second possibility we will describe an orientation 
of the emission plane by the azimuthal- and polar angle 
of the vector normal to that plane [27|. These angles 
are shown in figure |21 as <p* N and Q* N , respectively. Fur- 
ther the absolute orientation of the particles momenta 
in the emission plane will be described by ipN, the an- 
gle between rj meson and the vector product of the beam 
momentum and the vector N. 

Obviously, the interaction between the particles does 
not depend on the orientation of the emission plane, and 
therefore, it will fully manifest itself in the occupation 
density of the Dalitz plot which in our case will be repre- 
sented in terms of the square of the invariant masses of 
the two particle subsystems. The distribution of the ori- 
entation of the emission plane will reflect, however, the 
correlation between the intial and final channel and hence 
its determination should be helpful for the investigation 
of the production mechanism. 



III. EXPERIMENT AND DATA ANALYSIS 



Using the COSY-11 detection system |22j,|2j|, utilizing 
a stochastically cooled proton beam of the cooler syn - 
chrotron COSY [24| and a hydrogen cluster target [28| . 
we have performed a high statistics measurement of 
the pp — * ppr\ reaction at a nominal beam momentum 
of 2.027 GeV/c. The experiment was based on the 
four-momentum registration of both outgoing protons, 
whereas the r\ meson was identified via the missing mass 
technique. The positively charged particles have been 
identified combining the time of flight between the SI 
and S3 scintillation detectors and the momentum recon- 
structed by tracking trajectories registered by means of 
the drift chambers back to the target. The detection 
setup is sketched in figure 

The accuracy of the missing mass reconstruction de- 
pends on two parameters: the spread of the beam mo- 
mentum as well as the precision of the measured mo- 
menta of the outgoing protons. The latter is predomi- 
nantly due to the geometrical spread of the beam. Since 
for the reconstruction of the momenta we assume that the 
reaction takes place in the middle of the target we can- 
not correct on an event-by-event basis for the momentary 
spread of the beam. However, we can rectify the smear- 
ing due to the shifts of the centre of the beam relative 
to the target as well as the average changes of the ab- 
solute beam momentum during the experiment. There- 
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FIG. 3: Schematic view of the COSY-11 detection setup |22j. 
The cluster target [2s| is located in front of the accelerator 
dipole magnet. Positively charged particles which leave the 
scattering chamber through the thin exit foil are detected in 
two drift chamber stacks Dl, D2 and in the scintillator ho- 
doscopes SI, S2 and S3. Scintillation detector S4 and the 
position sensitive silicon pad detector Si are used in coinci- 
dence with the SI counter for the registration of the elastically 
scattered protons. Elastic scattering is used for an absolute 
normalisation of the cross sections of the investigated reac- 
tions and for monitoring both the geometrical spread of the 
proton beam and the position at which beam crosses the tar- 
get [23^ . In the left- lower corner a schematic view of the 
interaction region is depicted. 



fore, after the selection of events with two registered pro- 
tons, as a first step of the more refined analysis the data 
were corrected for the mean beam momentum changes 
(sec fig. determined from the measured Schottky fre- 
quency spectra and the known beam optics. In a next 
step, from the distributions of the elastically scattered 
protons, the Schottky frequency spectrum, and the miss- 
ing mass distribution of the pp — > ppX reaction, we have 
estimated that the spread of the beam momentum, and 
the spread of the reaction points in horizontal and verti- 
cal direction amount to a(pb eam ) = 0.63 ± 0.03 MeV/c, 
a(x) = 0.22 ± 0.02 cm, and <j{y) = 0.38 ± 0.04 cm, 
respectively. Details of this procedure can be found in ref- 
erences |23ll29j . Further on, a comparison of the experi- 
mentally determined momentum spectra of the elastically 
scattered protons with the distributions simulated with 
different beam and target conditions allows us to estab- 
lished the position at which the centre of the beam crosses 
the target with an accuracy of 0.25 mm [2]|. Accounting 
for the movement of the beam relative to the target we 
improved the missing mass resolution as demonstrated in 
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FIG. 4: The dashed curve denotes the proton beam mo- 
mentum distribution integrated over the whole measurement 
period. On the horizontal scale, the value of zero is set at a 
nominal beam momentum equal to 2.027 GeV/c. The solid 
line shows the beam momentum distribution after the correc- 
tion for the mean value which was determined in 10 second 
intervals. 
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FIG. 5: Missing mass distribution for the pp — > ppX reaction 
determined by means of the COSY-11 detection system at a 
beam momentum of 2.027 GeV. The solid histogram presents 
the data corrected for effects of the time dependent relative 
shifts between the beam and the target using the method 
described in reference [23j| . The dashed histogram shows the 
result before the correction. 



figure|3] After this betterment, the peak originating from 
the pp — > pprj reaction became more symmetric and the 
signal to background ratio increased significantly. The 
background was simulated taking into account pp — > ppX 
reactions with X — 2tt, Sir and 4ir. Since we consider 
here only the very edge of the phase space distribution 
where the protons are produced predominantly in the 
S-wave the shape of the background can be reproduced 
assuming that the homogenous phase space distribution 
is modified only by the interaction between protons. In- 
deed, as can be observed in figure HO the simulation de- 
scribes the data very well. The calculated spectrum is 
hardly distinguishable from the experimental points. 

The position of the peak on the missing mass spec- 
trum and the known mass of the r\ meson |30j enabled 
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FIG. 6: Missing mass spectrum for the pp — > ppX re- 
action determined in the experiment at a beam momen- 
tum of 2.0259 GeV/c. The mass resolution amounts to 
1 MeV/c 2 (<r). The superimposed histograms present the sim- 
ulation for 1.5 ■ 10 s events of the pp — > ppr) reaction, and 
10 10 events for the reactions pp — > pp2ir, pp — > pp3n and 
pp — > pp4ir. The simulated histograms were fitted to the data 
varying only the magnitude. The fit resulted in 24009 ± 210 
events with the production of the r\ meson. 



to determine the actual absolute beam momentum to 
be p beam = 2.0259 GeV/c ± 0.0013 GeV/c, which 
agrees within error limits with the nominal value of 
Pb°Z£ nal = 2 - 027 GeV/c. The real beam momentum 
corresponds to the excess energy of the pp?y system equal 
to Q = 15.5 ± 0.4 MeV. 



A. Covariance matrix and kinematical fitting 

As already mentioned in the previous section at the 
COSY-11 facility the identification of the pp — > pprj re- 
action is based on the measurement of the momentum 
vectors of the outgoing protons and the utilisation of 
the missing mass technique. Inaccuracy of the momen- 
tum determination manifests itself in the population of 
kinematically forbidden regions of the phase space, pre- 
venting a precise comparison of the theoretically derived 
and experimentally determined differential cross sections. 
Figure [7| visualizes this effect and clearly demonstrates 
that the data scatter significantly outside the kinemati- 
cally allowed region (solid line), in spite of the fact that 
the precision of the fractional momentum determination 
in the laboratory system {<j{pi a b)/piab ~ 7- 10~ 3 ) is quite 
high. Therefore, when seeking for small effects like for 
example the influence of the proton-ry interaction on the 
population density of the phase-space, one needs either 
to fold theoretical calculations with the experimental res- 
olution, or to perform the kinematical fitting of the data. 
Both procedures require the knowledge of the covariance 
matrix, and thus its determination constitutes a neces- 
sary step in the differential analysis and interpretation of 
the data. 

In order to derive the covariance matrix we need to 
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FIG. 7: Dalitz plot distribution of the pp — > ppr\ reaction 
simulated at Q = 15.5 MeV. The number of entries is shown 
in a logarithmic scale. The solid line gives the kinematically 
allowed area. The result was obtained taking into account the 
experimental conditions as described in the text. 



recognize and quantify all possible sources of errors in 
the reconstruction of the two proton momenta pi and p 2 . 
The four dominant effects are: i) finite distributions of 
the beam momentum and of the reaction points, ii) mul- 
tiple scattering in the dipole chamber exit foil, air, and 
detectors, iii) finite resolution of the position determina- 
tion of the drift chambers, and iv) a possible inadequate 
assignment of hits to the particle tracks in drift cham- 
bers in the case of very close tracks. Some of these, like 
the multiple scattering, depend on the outgoing protons' 
momenta, others, like the beam momentum distribution, 
depend on the specific run conditions and therefore must 
be determined for each run separately. 

In order to estimate the variances and covariances 
for all possible combinations of the momentum com- 
ponents of two registered protons we have generated 
1.5 • 10 s pp — > pprj events and simulated the response 
of the COSY-11 detection setup taking into account 
the above listed factors and the known resolutions of 
the detector components. Next, we analysed the sig- 
nals by means of the same reconstruction procedure as 
used in case of the experimental data. Covariances be- 
tween the i th and the j th components of the event vector 
(P = \pix,Piy,Piz,P2x,P2y,P2z]) were established as the 
average of the product of the deviations between the re- 
construced and generated values. The explicit formula 
for the sample of N reconstructed events reads: 
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FIG. 8: Variances of the protons's momentum components 
<y 2 (px), a 2 (p z ), and <J 2 (j> y ) shown as a function of the absolute 
values of the measured momenta. 



of the registered protons 2 there are only 12 independent 
values which determine the 6x6 error matrix V unam- 
biguously. 

Since inaccuracies of the momentum determination 
depend on the particle momentum itself (eg. multiple 
scattering) and on the relative momentum between pro- 
tons (eg. trajectories reconstruction from signals in drift 
chambers), we have determined the covariance matrices 
as a function of the absolute momentum of both protons: 
cov{i,j, \fi\, \p 2 \). 

As an example we present the covariance matrix for the 
mean values of \p~[\ and \p^\ in units of MeV 2 /c 2 , as es- 
tablished in the laboratory system with the z-coordinate 
parallel to the beam axis and y-coordinate corresponding 
to the vertical direction. 

P2x P2y P2z 

0.1 -3.0 
- -0.2 -0.2 
V = - - 37.0 - - 5.4 



j.reconJ ' 



fe=l 



Plx Ply Plz 

5.6 0.0 -13.7 1.7 

- 7.1 0.1 

- - 37.0 



Plx 
Ply 
Plz 
P2x 
P2y 
P2z 



(2) 



where P^ gen and P^ recon denote the generated and re- 
constructed values for the i th component of the vector P 
describing the k th event. 

Because of the inherent symmetries of the covariance 
matrix (cov(i,j) = cov(J,i)) and the indistinguishability 



2 The symmetry of all observables under the exchange of the two 
protons (pi <-> P2) implies that cov(i,j) = cov(i±3,j±3), where 
the '+' has to be taken for ij = 1,2,3 and the '-' for ij = 4,5,6. 
Thus for example cov(2,4) = cov(5,l). 
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FIG. 9: Spectrum of differences between generated and - af- 
ter simulation of the detector response - reconstructed abso- 
lute momenta of protons, as determined before (dashed line) 
and after the kinematical fit (solid line). 

The picture shows results obtained taking into account the 
experimental conditions as described in the text. 



Since the measurements have been performed close to 
the kinematical threshold the ejectile momentum com- 
ponent parallel to the beam is by far the largest one 
and its variance (var(p z ) = 37 MeV 2 /c 2 ) determines 
in first order the error of the momentum measurement. 
The second largest contribution stems from an anti- 
correlation between the x— and z— momentum compo- 
nents (cov(p x ,p z ) = — 13.7 MeV 2 /c 2 ), which is due 
to the bending of the proton trajectory - mainly in the 
horizontal direction - inside the COSY-11 dipole mag- 
net (see fig. |3J). There is also a significant correlation 
between the z components of different protons which is 
due to the smearing of the reaction points, namely, if 
in the analysis the assumed reaction point differs from 
the actual one, a mistake made in the reconstruction af- 
fects both protons similarly. Figure |HI depicts the varia- 
tion of var(p x ), var(p y ), and var{p z ) over the momen- 
tum plane (|pi|, \P2\)- Taking into account components 
of the covariance matrices V(|p!|, \p2\) and the distribu- 
tion of the proton momenta for the pp — > ppr/ reaction at 
Q = 15.5 MeV results in an average error for the measure- 
ment of the proton momentum of about 6 MeV/c. This 
can be also deduced from the distribution of the differ- 
ence between the generated and reconstructed absolute 
momenta of the protons. The corresponding spectrum is 
plotted as a dashed line in figure 

In the experiment we have measured 6 variables and 
once we assume that the event corresponds to the pp — > 
ppr\ reaction only 5 of them are independent. Thus we 
have varied the values of the event components demand- 
ing that the missing mass is equal to the mass of the 
r\ meson and we have chosen that vector which was the 
closest to the experimental one. The inverse of the co- 
variance matrix was used as a metric for the distance 



calculation. The kinematical fit improves the resolution 
by a factor of about 1.5 as can be concluded from com- 
paring the dashed- and solid lines in figureEl The finally 
resulting error of the momentum determination amounts 
to 4 MeV/c. 



B. Multidimensional acceptance corrections and 
results 



At the excess energy of Q = 15.5 MeV the COSY-11 
detection system does not cover the full 47r solid angle 
in the centre-of-mass system of the pp — > ppr\ reaction. 
Therefore, the detailed study of differential cross sections 
requires corrections for the acceptance. Generally, the 
acceptance should be expressed as a function of the full 
set of mutually orthogonal variables which describe the 
studied reaction unambiguously. As introduced in sec- 
tion [HI to define the relative movement of the particles 
in the reaction plane we have chosen two squares of the 
invariant masses: s pp and s p „, and to define the orien- 
tation of this plane in the center-of-mass frame we have 
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FIG. 10: Acceptance of the COSY-11 detection system for 
the pp — > ppr\ reaction at an excess energy of Q = 15.5 MeV 
presented as a function of s vv ,s VTl , and \cos(8^)\. The numbers 
were assigned to the bins in the three dimensional space s pp — 
s P n — \aos{6^)\ by first incrementing the index of s pp next of 
s pv and on the end that of |cos(f?^)|. Partitioning of | cos (6^)1 
into ten bins is easily recognizable. The two lower pictures 
show the acceptance for the first and the last bin of \cos(8^)\. 
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FIG. 11: Distribution of the polar angle of the emission of the 
r\ meson in the centre-of-mass system. Experimental data 
were corrected for the acceptance in the three dimensional 
space {s P p, s pv , \cos(0*)\). Full circles show the result with the 
assumption that the distribution of the ip angle is isotropic, 
and the open circles are extracted under the assumption that 
^ is as derived from the data (see text). Both results have 
been normalized to each other in magnitude. 
Open squares show (in arbitrary units) the distribution of 
the multi-pion production for the invariant mass of the pions 
ranging between ± 1 MeV around the mass of the r\ meson. 



taken the three Euler angles: The first two are simply 
the polar <p* and azimuthal 9* angles of the momentum 
of the r\ meson and the third angle ip describes the rota- 
tion of the reaction plane around the axis defined by the 
momentum vector of the rj meson. In the data evalua- 
tion we considerably benefit from the basic geometrical 
symmetries satisfied by the pp — > ppr/ reaction. Due to 
the axial symmetry of the initial channel of the two un- 
polarized colliding protons the event distribution over cj>* 
must be isotropic. Thus, we can safely integrate over 
</>*, ignoring that variable in the analysis. Furthermore, 
taking advantage of the symmetry due to the two iden- 
tical particles in the initial channel, without losing the 
generality, we can express the acceptance as a function 
of Spp,Spri,\cos(6*)\, and ip. To facilitate the calculations 
we have divided the range of \cos(9*)\ and ip into 10 bins 
and both s pp and s pri into 40 bins each. In the case 
of the Spp and s prj the choice was made such that the 
width of the interval corresponds to the standard devi- 
ation of the experimental accuracy. For |cos(#*)| and ip 
we have taken only ten partitions since from the previous 
experiments we expect only a small variation of the cross 
section over these variables [2(3, |2(| ■ In this represen- 
tation, however, the COSY- 11 detection system covers 
only 50% of the phase-space for the pp — > ppr/ reaction 
at Q = 15.5 MeV. To proceed with the analysis we as- 
sumed that the distribution over the angle ip is isotropic 
as it was for example experimentally determined for the 
pp — > ppuj, pp — > ppp, or pp — > pp<j) reactions [2^, l2t| . 
Please note that this is the only assumption of the reac- 
tion dynamics performed in the present evaluation. The 
validity of this supposition in the case of the pp — > ppr\ 
reaction will be discussed later. 
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FIG. 12: COSY-11 detection acceptance as a function of s pp 
and s pv , calculated under the assumption that the differential 
cross sections dc ^ r / g »- ) and ^ are isotropic. 



In the calculations we exploit also the symme- 
try of the cross sections under the exchange of the 
two identical particles in the final state which reads: 
vispirjjip) — &(s P2 r], ip + k). The resultant acceptance 
is shown in figure ITU1 One recognizes that only a small 
part (3%) of the phase-space is not covered by the de- 
tection system. In further calculations these holes were 
corrected according to the assumption of a homogeneous 
phase space distribution. Additionally it was checked 
that the corrections under other suppositions eg. regard- 
ing also the proton-proton FSI leads to negligible differ- 
ences. Full points in figure ^2 present the distribution 
of the polar angle of the r\ meson as derived from the 
data after the acceptance correction. Within the statis- 
tical accuracy it is isotropic. Taking into account this 
angular distribution of the cross section we can calculate 
the acceptance as a function of s pp and s pv only. This is 
shown in figure IT^l where one sees that now the full phase 
space is covered. This allows us to determine the distri- 
butions of Spp and s pri . The correctness of the performed 
procedures for the simulation of the detectors response, 
the event reconstruction programs, the kinematical fit- 
ting and acceptance correction can be confirmed by com- 
paring the distribution generated (figure 113b ) with the 
ones which underwent the complete analysis chain de- 
scribed in this section (figure EH 3 )- 

Knowing the distribution of the polar angle of the rj 
meson 9* and those for the invariant masses s pp and s pri 
we can check whether the assumption of the isotropy of 
the cross section distribution versus the third Euler's an- 
gle ip is corroborated by the data. For that purpose we 
calculated the acceptance as a function of ip and s prj as- 
suming the shape of the differential cross sections of 4 s — 
and (fcog^g*) a s determined experimentally. The obtained 

^ distribution is shown in figure IT4l and is not isotropic 
as assumed at the beginning. A fit of the function of 
the form ^ = a + b ■ \sin(tp)\ gives the value of 
b = 0.079 ± 0.014 fib/sr, indeed significantly different 
from the isotropic solution. This deviation cannot be 
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FIG. 13: (a) Dalitz plot distribution simulated for the 
pp —* pprj reaction at Q = 15.5 MeV. In the calculations 
the interaction between protons was taken into account, 
(b) Dalitz plot distribution reconstructed from the COSY-11 
detector response simulated for events from figure (a) taking 
into account the smearing of the beam and target, multiple 
scattering in the materials, and the detectors resolution. The 
evaluation included momentum reconstruction, kinematical 
fitting and the acceptance correction exactly in the same way 
as performed for the experimental data. The lines surround- 
ing the Dalitz plots depict the kinematical limits. 



assigned to any unknown behaviour of the background 
since the obtained distribution can be regarded as back- 
ground free. This is because the number of pp — > ppr\ 
events was elaborated for each invariant mass interval 
separately. The exemplarily missing mass spectra for the 
first, fourth, seventh and tenth interval of tp values, cor- 
rected for the acceptance are presented in figure IT5l 

From this figure one can infer that the shape of the 
background is well reproduced not only for the overall 
missing mass spectrum as shown previously in figure 
but also locally in each region of the phase space. Since 
the experimental data are quite well described by the sim- 
ulations we can rather exclude the possibility of a signif- 
icant systematical error which could cause the observed 
anisotropy of the differential cross section 9S. ■ 

The evaluated distribution is however in disagreement 
with our working assumption that the cr(tp) is isotropic. 
Therefore we performed a full acceptance correction pro- 
cedure from the very beginning assuming that the dis- 
tribution of 4t is as determined from the data. After 
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FIG. 14: Distribution of the cross section as a function of 
the angle ip as determined in the first iteration. The su- 
perimposed histogram corresponds to the fit of the function 
jfe = a + b ■ \sin(ip)\. The range of the %ji angle is shown 
from to 7T only, since in the analysis we take advantage of 
the symmetry = + 
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FIG. 15: Missing mass distributions for the first, fourth, sev- 
enth and tenth bin of tp with the superimposed lines from 
the simulation of pp — » pprj and the multi-pion background 
pp — » ppljrm) reactions. Amplitudes of simulated distribu- 
tions were fitted to the experimental points. 



repeating the procedure three times we observed that 
the input and resultant distributions are in good agree- 
ment. The result after the third iteration is shown in 
figure [TBI bv the full circles. It is only slightly different 
from the one obtained after the first iteration as shown 
in figure To rise the confidence of the convergence of 
the performed iteration we accomplished the full proce- 
dure once more, but now assuming that the distribution 
of ^ is much more anisotropic than determined from the 
data. As an entry distribution we took the dashed line 
shown in figure 1161 Again after two iterations we have 
got the same result as before. Finally determined values 
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FIG. 16: Distribution of the cross section as a function 
of the angle ip. Full circles stand for the final results 
of the obtained after three iterations. The superim- 
posed histogram (solid line) corresponds to the fit of the 
function ^ = a + b ■ \sin(ip)\ which resulted in 
a = 0.186 ± 0.004 fib/sr and b = 0.110 ± 0.014 fib/sr. 
The dashed line shows the entry distribution used for the sec- 
ond series of iterations as described in the text. 
Open circles represent the data from the left upper corner of 
the Dalitz plot (see for example figure IT^jl . At that region 
of the Dalitz plot due to the non-zero four dimensional ac- 
ceptance over (s pp , s pv , \cos(8*)\,ip) bins the spectrum (open 
circles) was corrected without a necessity of any assumptions 
concerning the reaction cross section. 



TABLE I: Differential cross section in the ip angle for the 
pp — » ppr\ reaction at Q = 15.5 MeV. Fitting to the data 
a function of the form & — a + b \sin(ip)\ resulted in 
a = 0.186 ± 0.004 fib/sr and b = 0.110 ± 0.014 fib/sr. 
Figure IHTl illustrates the result. 



iP[rad] g 


fib 


da 


iib 




rad 


dip 


rad 





experiment fit 

0157 0.200 ± 0.004 O202 

0.471 0.249 ± 0.014 0.235 

0.785 0.269 ± 0.021 0.263 

1.100 0.242 ± 0.025 0.283 

1.414 0.309 ± 0.031 0.294 

1.728 0.223 ± 0.025 0.294 

2.042 0.334 ± 0.031 0.283 

2.356 0.268 ± 0.023 0.263 

2.670 0.281 ± 0.015 0.235 

2.985 0.204 ± 0.004 0.202 



for ^("0) are given in tabled To corroborate this obser- 
vation we have evaluated the distribution over ip angle 
(see figure ITB|) from the phase space region which has no 
holes in the acceptance expressed as a four-dimensional 
function of the variables s pp , s pv , \cos(9*)\, and ip, this is 
for the values of s pp and s pv corresponding to the upper 
left corner of figure ^] Again the obtained distribution 
presented as open circles in figure IT^l is anisotropic, and 
moreover agrees with the spectrum determined from all 
events. The anisotropy of the cross section in the ip an- 



gle reflects itself in an anisotropy of the orientation of the 
emission plane. 

The determined cross section distribution in function 
of the polar angle 9* N of the vector normal to that plane 
is shown in figure 1171 and the corresponding values are 
listed in table[n] The distribution is not isotropic, which 
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FIG. 17: Differential cross section as a function of the polar 
angle of the vector normal to the emission plane. 

is particularly visible for the low values of |cos(#^)| bur- 
dened with small errors. As depicted in figure the 
\cos(9* N )\ = denotes such configuration of the ejectilcs 
momenta in which the emission plane comprises the beam 
axis. In that case the acceptance of the COSY-11 de- 
tection system is much larger than for the configuration 
where the emission plane is perpendicular to the beam. 
Due to this reason the error bars in figure 1171 increase 
with growth of \cos(9^)\. It is worth to stress that the 
tendency of the ppr\ system to be produced preferentially 
if the emission plane is perpendicular to the beam is in 
line with the preliminary analysis of the experiment per- 
formed by the TOF collaboration |23|. Elucidation of 
that non-trivial behaviour can reveal interesting features 
of the dynamics of the production process. 

It is important to note that the shape of the s pp , s mi 
and cos(9*) distributions keeps unchanged during the 
whole iteration procedure. Figure EH shows the distri- 
butions of the square of the proton-proton and proton-?? 
invariant masses. The spectra after the second and third 
iterations are shown. One recognizes that the form of the 
spectra remains unaltered. The same conclusion can be 
drawn for the da , a ,^> distribution as demonstrated in 
figure ITTI From that comparison one can conclude that 
the shapes of the determined distributions are - within 
the statistical accuracy - independent of the shape of the 
cross section and can be treated as derived in a com- 
pletely model independent manner. Similarly, as in the 
case of the ^ distribution, the differential cross sections 
in all other variables reported in this article are not dete- 
riorated by the background. This is because the number 
of pp — y ppr] events was determined for each investigated 
phase-space interval separately. As an example the miss- 
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FIG. 18: Distributions of the invariant masses s pp and s pv 
determined for the two different assumptions about the cross 
section dependence of the V angle. 



ing mass spectra for three bins of the proton-?? invariant 
mass are presented in figure 1191 As already noticed for 
^ the shape of the background is well reproduced locally 
in each region of the phase space. 
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FIG. 19: Missing mass distributions determined for the in- 
variant mass bins as depicted inside the figures. The spectra 
were corrected for the acceptance. The histograms show the 
simulations for the multi-pion background pp — > pp(m-K) and 
the pp — > pprj reactions fitted to the data with the amplitudes 
as the only free parameters. 



Total and differential cross sections 



Though the form of the s pp , s PV , and cos(8*) is in- 
dependent of the ip distribution, the total cross section 
derived from the data depends on the shape of ^ quite 
significantly. It amounts to 3.24 ± 0.03 fib, yet it changes 
by ± 0.2 fib when varying the parameters of the function 
^ = a + b ■ \sin(ip)\ by ± three standard deviations. 
Therefore, we use that variation as an estimation of the 
systematic error in the acceptance correction. To this we 
must add a 3% systematical uncertainty stemming from 
the luminosity determination pflj . The luminosity was 
determined by the comparison of the measured differen- 
tial distribution of the elastically scattered protons with 
the results of the EDDA collaboration [3l|]. The deter- 
mined value amounts to 811 ± 8 ± (3%) nb^ 1 . Thus 
the overall systematical error of the cross section value 



amounts to 0.30 fib. In summary, we determined that at 
Q = 15.5 MeV the total cross section for the pp — > pprf 
reaction is equal to 3.24 ± 0.03 ± 0.30 fib, where the first 
and second error denote the statistical and systematical 
uncertainty, respectively |3^ . 

In the below tables and figures we will present the val- 
ues for the differential cross section concerning the ob- 
servables described in section [n] If possible the data will 
be compared to the result of measurements performed at 
the non-magnetic spectrometer COSY-TOF [2£j ■ An in- 
terpretation of the elaborated distribution follows in the 
next section. The distributions of the squared invariant 
masses as presented already in figure ^1 are listed in ta- 
blc lIIII The distribution of the polar angle of the rf meson 
emission in the centre-of-mass system is given in table [H] 
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TABLE II: Differential cross section in \cos(8^)\, \cos(9^)\, 
\cos(8p P )\, and \cos(8pp)\ for the pp — > ppr\ reaction at 
Q = 15.5 MeV. 
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FIG. 20: Differential cross section of the pp — > ppr\ reaction 
as a function of the rj meson centre-of-mass polar angle. 
Full circles depict experimental results for the pp —* ppr\ re- 
action measured at Q = 15.5 MeV by the COSY-11 collabo- 
ration (this article) and the open circles were determined by 
the TOF collaboration at Q = 15 MeV Hj. The TOF points 
were normalized in amplitude to our result, since for that data 
the absolute scale is not evaluated. 



and shown in figure |2*U1 Clearly, our data agree very well 
with the angular dependence determined by the TOF col- 
laboration. As already mentioned, the two identical par- 
ticles in the initial state of the pp — > ppr\ reaction imply 
that the angular distribution of either ejectile must be 
symmetric around 90 degree in the centre-of-mass frame. 
Here we use that reaction characteristic, yet in the pre- 
vious analysis of this data [3i| we utilized this property 
to check the correctness of the acceptance calculation. In 
reference [33] we presented the differential cross section of 
the 77 meson centre-of-mass polar angle for the full range 
of cos(9*) and found that this is completely symmetric 
around cos(9*) — 0. 

For completeness we calculated also a distribution of 
the angle ipN defining the orientation of the ppr\ system 
within the emission plane. The result is shown in fig- 
ure [2H and the values of -4f- are listed in table HVl 

Up to now we presented invariant mass spectra of the 



proton-proton and proton-77 systems and angular distri- 
bution for two sets of non-trivial angles which describe 
the orientation of ejectiles within the emission plane and 
the alignment of the plane itself, namely (cos(#*), ip) and 
(cos{9* N ), ip N ). 

Since one of the important issues which we will dis- 
cuss in the next section is the contribution from higher 
partial waves we evaluated also an angular distribution 
of the relative momentum of two protons seen from the 
proton-proton centre-of-mass subsystem (see figure 1221) . 
The distribution of that angle should deliver information 



TABLE III: Distribution of the square of the invariant mass of 
the proton-proton (s pp ) and proton-?; {s pv ) systems measured 
at Q = 15.5 MeV via pp — > ppr\ reaction. 



[GeV 2 /c 



3.5221 
3.5236 
3.5251 
3.5265 
3.5280 
3.5294 
3.5309 
3.5324 
3.5338 
3.5353 



3.5367 
3.5382 
3.5397 
3.5411 
3.5426 
3.5440 
3.5455 
3.5469 
3.5484 
3.5499 



3.5513 
3.5528 
3.5542 
3.5557 
3.5572 
3.5586 
3.5601 
3.5615 
3.5630 
3.5645 



3.5659 
3.5674 
3.5688 
3.5703 
3.5718 
3.5732 
3.5747 
3.5761 
3.5776 
3.5790 



da 

UiSpp 

fib 



GeU 2 /c 4 



48.7 ± 1.7 

95.8 ± 3.0 
111.1 ± 3.6 



96.7 
94.1 
84.9 
76.7 
79.4 
71.3 
70.4 



3.8 
4.4 
4.0 
4.0 
4.3 
4.3 
4.4 



77.4 ± 4.7 



62.5 
58.7 
66.7 
57.8 
64.4 
65.0 
54.6 
61.3 
58.4 



4.5 
4.5 
4.8 
4.9 
5.1 
5.1 
4.7 
5.1 
5.5 



52.4 ± 5.2 

50.7 ± 5.1 

46.8 ± 5.1 

51.5 ± 5.4 

46.8 ± 5.3 
50.7 ± 5.5 

50.9 ± 5.5 
48.3 ± 5.7 
45.7 ± 5.7 
39.3 ± 5.2 



39.0 ± 5.2 



± 
± 
± 
± 
± 
± 

26.0 ± 3.4 
16.3 ± 1.3 
7.7 ± 1.4 



33.6 
33.0 
34.4 
31.8 
27.9 
31.9 



4.7 
4.6 
4.6 
4.4 
4.2 
3.9 



[GW 2 /c 4 ] 



2.2075 
2.2087 
2.2098 
2.2110 
2.2121 
2.2133 
2.2144 
2.2156 
2.2168 
2.2179 



2.2191 
2.2202 
2.2214 
2.2225 
2.2237 
2.2248 
2.2260 
2.2272 
2.2283 
2.2295 



2.2306 
2.2318 
2.2329 
2.2341 
2.2353 
2.2364 
2.2376 
2.2387 
2.2399 
2.2410 



2.2422 
2.2433 
2.2445 
2.2457 
2.2468 
2.2480 
2.2491 
2.2503 
2.2514 
2.2526 



da 

dSp-q 



GeV 2 /c 



16.8 
27.9 
26.2 
33.5 
36.2 
44.3 
43.6 
51.6 
56.1 
50.2 



1.4 
1.7 

2.3 
2.8 
3.1 
3.3 
3.5 
3.8 
3.9 
4.0 



57.6 
56.9 
67.9 
65.8 
69.8 
70.3 
81.7 
68.0 
72.9 
89.0 



4.2 
4.2 
4.2 
4.2 
4.4 
4.4 
4.7 
4.2 
4.3 
4.7 



80.3 ± 4.5 

84.5 ± 4.5 

83.6 ± 4.4 
95.6 ± 4.6 
100.2 ± 4.7 

93.6 ± 4.4 

99.7 ± 4.5 
102.6 ± 4.3 
103.6 ± 4.3 
104.5 ± 4.1 



91.9 
96.1 
97.8 
94.5 
92.3 
88.1 
73.3 
65.7 
45.7 
22.2 
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1.4 
0.7 
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FIG. 21: Differential cross section in t/>jv for the pp —* ppr\ re- 
action measured at Q = 15.5 MeV. The variable i/jn is defined 
in section ITTl 



about the partial waves of the proton-proton system in 
the exit channel. In case of the two body scattering, the 
beam line, which is at the same time the line along which 
the centre-of-mass system is moving, constitutes the ref- 
erence frame for the angular distributions. For a three 
body final state the beam axis is not a good direction 
to look for the angular distributions relevant for the 



TABLE IV: Differential cross section in ipN for the pp — > 
pprj reaction measured at Q = 15.5 MeV. The variable ip^ is 
defined in section [2] Values are presented in figure |2*T1 



ipN [rad] 


da 

dip N 


fib 

sr 


0.314 


0.242 ± 0.007 


0.942 


0.261 ± 0.008 


1.571 


0.277 ± 0.007 


2.199 


0.261 ± 0.008 


2.827 


0.255 ± 0.007 


3.456 


0.267 ± 0.009 


4.084 


0.254 ± 0.009 


4.712 


0.255 ± 0.009 


5.341 


0.251 ± 0.009 


5.969 


0.263 ± 0.008 



relative angular momenta of the two particles |34] . But 
by analogy to the two body system, an instructive refer- 
ence axis for angular distributions in the proton-proton 
subsystem is now the momentum of the recoil rf meson, 
since the direction of that meson is identical with the 
direction of the movement of the proton-proton centre- 
of-mass subsystem. The distribution of the differential 
cross section in cos(9**) is given in table ITTl In this table 
we listed also the differential cross section as a function 
of angle 9* p of the relative proton momentum seen from 
the overall centre-of-mass frame, as this is often consid- 
ered in the theoretical works. In figure 1251 our results are 
compared to the angular distribution extracted by the 
TOF collaboration. Both experiments agree very well 




FIG. 22: Definition of 8pp, the polar angle of the relative 
proton-proton momentum with respect to the momentum of 
the rj meson as seen in the di-proton rest frame. Lower picture 
shows differential cross section in #** as determined for the 
pp — > pprf reaction at Q = 15.5 MeV. 



within the statistical accuracy, and both indicate a slight 
decrease of the cross section with increasing \cos(9* p )\. 



IV. INTERPRETATION OF RESULTS 

The interaction between particles depends on their rel- 
ative momenta or equivalently on the invariant masses of 
the two-particle subsystems. Therefore it should show 
up as modification of the phase-space abundance in the 
kinematical regions where the outgoing particles possess 
small relative velocities. Only two invariant masses of 
the three subsystems are independent and therefore the 
entire accessible information about the final state interac- 
tion of the three-particle system can be presented in the 
form of the Dalitz plot. The upper panel of figurel24lindi- 
cates the event distribution as determined experimentally 
for the pprf system at an excess energy of Q = 15.5 MeV. 
In this figure one easily recognizes the growths of the 
population density at the region where the protons have 
small relative momenta which can be assigned to the 
strong attractive S-wave interaction between the two pro- 
tons. This is qualitatively in agreement with the ex- 
pectation presented in figure I24f middle), which shows 
the results of Monte-Carlo calculations where the homo- 
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FIG. 23: Distribution of the centre-of-mass polar angle of the 
relative protons momentum with respect to the beam direc- 
tion determined for the pp — > ppr\ reaction at Q = 15.5 MeV. 
The COSY-11 result (closed circles) is compared to the data 
points determined at Q = 15 MeV by the TOF collaboration 
(open circles) |2Cj . 



geneously populated phase space was weighted by the 
square of the on-shell 1 Sq proton-proton scattering am- 
plitude. However, already in this two dimensional repre- 
sentation it is visible that the experimentally determined 
distribution remains rather homogeneous outside the re- 
gion of the small proton-proton invariant masses, whereas 
the simulated abundance decreases gradually with grow- 
ing s pp (as indicated by the arrow). The lower panel of 
figure |^] shows the simulated phase-space density dis- 
tribution disregarding the proton-proton interaction but 
accounting for the interaction between the ?7-meson and 
the proton. Due to the lower strength of this interac- 
tion the expected deviations from a uniform distribution 
is smaller by about two orders of magnitude, yet an en- 
hancement of the density in the range of low invariant 
masses of proton-ry subsystems is clearly visible. Note 
that the scale in the lower figure is linear whereas in the 
middle and upper panel it is logarithmic. Due to weak 
variations of the proton-?7 scattering amplitude the en- 
hancement originating from the ?7-meson interaction with 
one proton is not separated from the ?7-meson interac- 
tion with the second proton. Therefore an overlapping of 
broad structures occurs. It is observed that the occupa- 
tion density grows slowly with increasing s pp , opposite to 
the effects caused by the S-wave proton-proton interac- 
tion, but similar to the modifications expected for the NN 
P-wave |35| . From the above example it is obvious that, 
only from experiments with high statistics, signals of the 
meson-nucleon interaction can be observed over the over- 
whelming nucleon-nucleon hnal state interaction. A de- 
viation of the experimentally observed population of the 
phase-space from the expectation based on the mentioned 
assumptions is even better visible in figure E5l This figure 
presents the projection of the phase-space distribution 
onto the s pp axis corresponding to the axis indicated by 
the arrows in the two lower parts of ngure Ell The super- 
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FIG. 24: Dalitz plot distributions. 

(upper panel) Experimental result determined for the pp — > 
pprj reaction at Q = 15.5 MeV. Data were corrected for the 
detection acceptance and efficiency. 

(middle) Monte-Carlo simulations for the pp —* pprj reaction 

at Q = 15.5 MeV: Phase-space density distribution modified 

by the 1 5'o proton - proton final state interaction. 

(lower panel) Simulated phase-space density distribution 

modified by the proton-77 interaction with a scattering length 

of a pv — 0.7 fm + i 0.3 fm. Details of the calculations together 

with the discussion of the nucleon-nucleon and nucleon-meson 

final state interaction can be found in reference fl4| . 

The scale in the figure is linear in contrast to the above panels. 

Lines surrounding the plots show the kinematical limits. 
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FIG. 25: Distributions of the square of the proton- 
proton (s pp ) and proton-?? (s pv ) invariant masses determined 
experimentally for the pp —> ppq reaction at the excess energy 
of Q = 15.5 MeV by the COSY-11 collaboration (closed cir- 
cles), at Q = 15 MeV by TOF collaboration (open circles) [2o|. 
and at Q = 16 MeV by PRO MICE/ WAS A (open trian- 
gles) 0. The TOF and PROMICE/WASA data have been 
normalized to those of COSY-11, since these measurements 
did not evaluate the luminosities but rather normalized the re- 
sults to reference [f| (see also comment |3^|). The integrals of 
the phase space weighted by the square of the proton-proton 
on-shell scattering amplitude (dotted lines)-FSI pp , and by the 
product of FSIpp and the square of the proton-77 scattering 
amplitude (thick solid lines), have been normalized arbitrar- 
ily at small values of s pp . The thick solid line was obtained 
assuming a scattering length of a pv — 0.7 fm + i 0.4 fm. 
The expectation under the assumption of the homogeneously 
populated phase space are shown as thin solid curves. 



imposed lines in figure I5B1 correspond to the calculations 
performed under the assumption that the production am- 
plitude can be factorized into a primary production and 
final state interaction. The dotted lines result from cal- 
culations where only the proton-proton FSI was taken 
into account, whereas the thick-solid lines present results 
where the overall enhancement was factorized into the 
corresponding pair interactions of the ppr) system. This 
factorisation Ansatz is of course only valid if the different 
amplitudes are completely decoupled which is certainly 
not the case here. Therefore, this calculations should be 
considered as a rough estimate of the effect introduced 
by the FSI in the different two body systems. The en- 
hancement factor accounting for the proton-proton FSI 
has been calculated [T3 . l38| as the square of the on-shell 
proton-proton scattering amplitude derived according to 
the modified Cini-Fubini-Stanghellini formula including 
the Wong-Noyes Coulomb corrections [3t|. The homo- 
geneous phase-space distribution (thin solid lines) devi- 
ate strongly from the experimentally determined spectra. 
The curves including the proton-proton and proton-77 FSI 
reflect the shape of the data for small invariant masses of 
the proton-proton system, yet they deviate significantly 
for large s pp and small s pri values. An explanation for 
this discrepancy could be a contribution from P-wave 
proton-proton interaction [Tlj . or a possibly inadequate 
assumption that proton-77 and proton-proton interaction 
modify the phase space occupations only as incoherent 
weights [H. 

Slightly better description is achieved when the proton- 
proton interaction is accounted for by the realistic 
nucleon-nucleon potential. The upper picture in figure I2tj| 
depicts the results obtained using two different models 
for the p roducti on p rocess as well as for the NN interac- 
tion (TflEU UlEa. The calculations for the 3 P ->%s 
transition differ slightly, but the differences between the 
models are, by far, smaller than the observed signal. 

Therefore we can safely claim that the discussed effect 
is rather too large to be caused by the particular assump- 
tions used for the production operator and NN potential. 

As can be seen in figure [501 the distribution of the 77 
polar angle in the center of mass frame is fully isotropic. 
This is the next evidence — besides the shape of the 
excitation function and the kinematical arguments dis- 
cussed in reference " ~ that at this excess energy 
(Q = 15.5 MeV) the 77 meson is produced in the center 
of mass frame predominantly with the angular momen- 
tum equal to zero. Similarly, the distribution determined 
for the polar angle of the relative proton-proton momen- 
tum with respect to the momentum of the 77 meson as 
seen in the di-proton rest frame is also consistent with 
isotropy. Anyhow, even the isotropic distribution in this 
angle does not imply directly that the relative angular 
momentum between protons is equal to zero, because of 
their internal spin equal to | . Therefore, the contribution 
from the 3 P -wave produced via the 1 Sq ^ 3 PqS transi- 
tion cannot be excluded. The isotropic angular distribu- 
tion, as pointed out in reference |ll| . can also be prin- 
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FIG. 26: (upper picture) Distribution of the square of the 
proton-proton (s pp ) invariant mass for the pp — > pprj reac- 
tion at an excess energy of Q = 15.5 MeV. Solid and dashed 
lines corresponds to the calculations under the assumption of 
the 3 Po — ^Sos transition according to the models described 
in references |TcL l41| and [TH E^. respectively. The dotted 
curve shows the result with the inclusion of the 1 So — > 3 Pos 
contribution as suggested in reference [Tl|. 
(lower picture) The same data as above but with curves 
denoting preliminary three-body calculations [l^| of the fi- 
nal pprj system as described in |4p|l . At present only the 
dominant transition 3 Po — Sos is taken into account and 
the production mechanism is reduced to the excitation of 
the 511(1535) resonance via the exchange of the 7r and r\ 
mesons. The solid line was determined with the rigor- 
ous three-body approach [l8| where the proton-proton sec- 
tor is described in terms of the separable Paris potential 
(PEST3) EJJ, and for the 77-nucleon scattering amplitude an 
isobar model analogous to the one of reference |44| is used 
with a V M = 0.5 fm + i ■ 0.32 fm. The dashed line is ob- 
tained if only pairwise interactions (pp + prj) are allowed. The 
effect of proton-proton FSI at small s pp is overestimated due 
to neglect of Coulomb repulsion between the protons. The 
lines are normalized arbitrarily but their relative amplitude 
is fixed from the model. 



TABLE V: Square of the proton-proton invariant mass from 
the pp — > pprj reaction measured for the excess energy range 
4 MeV < Q < 5 MeV HEJ. 



s pp [GeV jc 



da 
ds 



3.52164 
3.52204 
3.52264 
3.52344 
3.52444 
3.52564 
3.52704 
3.52864 
3.53044 
3.53244 
3.53464 
3.53704 
3.53964 



££_ 



fib 



GeV 2 /c 4 



19.0 
74.0 
56.9 
91.7 
69.8 
56.9 



13 
21 
15 
17 
13 
11 



51.2 ± 9.9 
41.8 ± 8.4 
43.4 ± 8.1 

36.6 ± 7.0 

24.7 ± 5.5 
6.8 ± 2.8 
1.1 ± 1.1 



cipally achieved by the destructive interference between 
the transitions 1 Sq —^ 3 PqS and 1 D 2 — > 3 i"2S- 

As show in reference [Tl"| the invariant mass distribu- 
tions can be very well described when including higher 
partial-wave amplitudes. In fact, as depicted by the dot- 
ted line in the upper panel of figure EEl an admixture of 
the 1 5o ^ 3 PqS transition leads to the excellent agree- 
ment with the experimentally determined invariant mass 
spectra. However, at the same time, the model of refer- 
ence Ej| leads to strong discrepancies in the shape of the 
excitation function as can be deduced from the compar- 
ison of the dashed-dotted line and the data in figure ^ 
Whereas it describes the data points in the excess energy 
range between 40 MeV and 100 MeV it underestimates 
the total cross section below 20 MeV by a factor of 2. 

Interestingly, the enhancement at large s pp is visible 
also at much lower excess energy. This can be con- 
cluded from figure 1771 in which the COSY-11 data at 
Q ps 4.5 MeV [46| are compared to the simulations based 
on the assumption that the phase-space abundance is due 
to the proton-proton FSI only. This observation could 
imply that the effect is caused by the proton-77 interac- 
tion rather than higher partial waves, since their contri- 
bution at such small energies is quite improbable [l4j| . 
However, as shown in the lower part of figure 1251 the rig- 
orous three-body treatment of the pprj system leads at 
large values of s pp to the reduction of the cross section in 
comparison to the calculation taking into account only 
first order rescattering (pp+pr?) 18]. Here, both calcu- 
lations include only the 3 Pq —* 1 SqS transition. Though 
the presented curves are still preliminary, we can quali- 
tatively assess that the rigorous three-body approach, in 
comparison to the present estimations, will on the one 
hand enhance the total cross section near threshold as 
shown in figure ^ while on the other hand it will de- 
crease the differential cross section at large values of s pp . 
This is just opposite to the influence of P-waves in the 
proton-proton system. 
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FIG. 27: Distribution of square of the proton-proton invari- 
ant mass from the pp — * pprj reaction measured at COSY-11 
for the excess energy range 4 MeV < Q < 5 MeV 0, . Nu- 
merical values are listed in table [3 The superimposed line 
shows the result of simulations performed under the assump- 
tion that the phase space population is determined exclusively 
by the on-shell interaction between outgoing protons. The 
"tail" at large s pp values is due to the smearing of the ex- 
cess energy, since this former COSY-11 data have not been 
kinematically fitted. Additionally to the 1 MeV range of Q a 
smearing of about 0.3 MeV (a) should be taken into account. 



From the above presented considerations it is rather 
obvious that the rigorous three-body treatment of the 
produced pprj system and the exact determination of the 
contributions from the higher partial waves may result in 
the simultaneous explanation of both observations: the 
near-threshold enhancement of the excitation function of 
the total cross section, and the strong increase of the 
invariant mass distribution at large values of s pp . For 
the unambiguous determination of the contributions from 
different partial waves spin dependent observables are re- 
quired [lj . The first attempt has been already reported 



V. CONCLUSION 

Using the stochastically cooled proton beam at the 
Cooler Synchrotron COSY and the COSY-11 facility 



we have determined the total and differential cross sec- 
tions for the pp — > pprj reaction at an excess energy of 
Q = 15.5 MeV. The high statistics data sample allowed 
us for the clear separation of events corresponding to the 
pp — » ppr\ reaction from the multi-pion production at 
each investigated phase-space bin, and the multidimen- 
sional acceptance correction allowed to extract the result 
without necessity of any assumption about the reaction 
process. 

The determined distributions of the centre-of-mass po- 
lar angle of the r\ meson emission as well as the distri- 
bution of the relative proton-proton momentum with re- 
spect to the momentum of the r\ meson are consistent 
with isotropy. Though, in the latter a small tendency of 
an increase of the cross section at 90 degree is observed. 
In contrary a rather strong decrease of the cross section 
was found at 90 degree for the center-of-mass polar angle 
of the vector normal to the emission plane. Explanation 
of that effect may reveal an interesting characteristic of 
the dynamics of the production process. 

The determined invariant mass spectra of the two- 
particle subsystems deviate strongly from the predictions 
based on the homogenous population of events over the 
phase space. Deviations at low proton-proton invariant 
mass values can be well explained as an influence of the 
S-wave interaction between the two protons. However, 
an unexpectedly large enhancement of the occupation 
density in the kinematical regions of low 77-proton rela- 
tive momentum is not yet understood. We have demon- 
strated that for the simultaneous description of the exci- 
tation function and invariant mass distributions a rigor- 
ous three-body calculation with inclusion of the contri- 
bution from higher partial waves is needed. 
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